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ABSTRACT 

We present axisymmetric maximum disk mass models for a sample of 74 spiral galaxies taken from 
the southern sky Fabry-Perot TuUy-Fisher survey (Schommer et al. 1993). The sample contains galaxies 
spanning a large range of morphologies and having rotation widths from 180 km s^^ to 680 km s^^. For 
each galaxy we have an J-band image and a two dimensional Ha velocity field. We decompose the disk 
and bulge by fitting models directly to the /-band image. This method utilizes both the distinct surface 
brightness profiles and shapes of the projected disk and bulge in the galaxy images. The luminosity 
profiles and rotation curves are derived using consistent centers, position angles, and inclinations derived 
from the photometry and velocity maps. The distribution of mass is modeled as a sum of disk and bulge 
components with distinct, constant mass-to-light ratios. No dark matter halo is included in the fits. 

The models reproduce the overall structure of the rotation curves in the majority of galaxies, providing 
good fits to galaxies which exhibit pronounced structural differences in their surface brightness profiles. 
75% of galaxies for which the rotation curve is measured to R23.5 or beyond are well fit by a mass-traces- 
light model for the entire region within R23.5. The models for about 20% of the galaxies do not fit well; 
the failure of most of these models is traced directly to non-axisymmetric structures, primarily bars but 
also strong spiral arms. The median /-band M/L of the disk plus bulge is 2.4 ± O.Ohys in solar units, 
consistent with normal stellar populations. These results require either that the mass of dark matter 
within the optical disk of spiral galaxies is small, or that its distribution is very precisely coupled to the 
distribution of luminous matter. 



1. INTRODUCTION 

Extended H I rotation curves provide deeply compelling 
evidence for a dark matter component that dominates the 
total mass of spiral galaxies. However, the fraction of 
lunnnous to dark matter (L/D) within the optical disk 
(within the 25 mag/arcsec^ i3-band, 23.5 mag/arcsec^ /- 
band isophote) is very poorly known. Typically, extended 
rotation curves have no distinct features which might re- 
flect the end of the luminous disk (Bahcall & Casertano 
1985). As a consequence, they provide virtually no con- 
straint on L/D within the optical disk. Mass models in 
which the dark component dominates the inner mass dis- 
tribution generally fit the rotation curves within the opti- 
cal disk as well as models with a negligible dark component 
(van Albada & Sancisi 1986, Lake & Feinswog 1989). 

The mass and luminosity of spiral galaxies are, however, 
very strongly connected. Spiral galaxies exhibit a very 
tight, one parameter relation between rotational velocity 
and luminosity, the Tully-Fisher (TF) relation (TuUy & 
Fisher 1977, see Jacoby et al. 1992 for a review). The to- 
tal mass-to-light ratio (M/L) within the optical disk is re- 
markably uniform among all spirals (Rubin 1985, Roberts 
& Haynes 1994) including low surface brightness spirals 
(Sprayberry et al. 1995). Indeed, the typical value of 
M/L (we derive here (M/L)/ = 2.4±0.9) is consistent with 



that expected from normal stellar populations (Larson & 
Tinsley 1978, Bruzual & Chariot 1993, Worthey 1994), 
with no dark matter. Optical rotation curves are not, in 
general, fiat and featureless. They span a range of shapes 
from linearly rising to falling with radius (Rubin 1985) and 
also display smaller scale "bumps and wiggles" (Freeman 
1992). Mass models fitted to rotation curves within the 
optical disk rarely require dark matter halos to yield good 
fits (Kalnajs 1983, Kent 1986, Buchhorn 1992). These 
models generally reproduce the large scale features and 
sometimes reproduce the smaller scale "bumps and wig- 
gles" in the optical rotation curves. 

Dynamical arguments suggest that the inner regions of 
spirals cannot be dominated by dark matter. Within the 
optical radius, a dark matter halo with an average pro- 
jected surface mass density greater than that of the disk 
would act to suppress common instabilities such as bars 
and 2-arm spiral structure (Athanassoula et al. 1987). 
The existence of lopsided modes, which appear in many 
disk galaxies (Rix & Zaritsky 1995), may require even 
lower halo mass densities. A bulge or thick disk also act to 
stabilize the disk, further lowering the allowed dark mat- 
ter density. It is possible to contsruct a model disk galaxy 
with an almost flat rotation curve which is stable with no 
dark matter (Sellwood & Evans 2000). This reverses the 
disk-stability argument Ostriker & Peebles (1973), which 
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is often used to support the need for a massive dark mat- 
ter halo within the optical radius. Of course, spherical 
halo mass outside the optical radius has no effect on disk 
stability. A stellar bar would also interact with a massive 
halo through dynamical friction. For large dark matter 
densities this interaction acts to rapidly slow the pattern 
rotation speed of the bar (Weinberg 1985, Debattista & 
Sellwood 2000) far below that expected in real galaxies 
(Sellwood & Wilkinson 1993), or measured in any galaxy 
(Merrifield & Kuijken 1995, Gerssen, Kuijken & Merrifield 
1999) 

Simulations of galaxy formation through dissipationless 
collapse of matter yield strongly triaxial halos (Dubinski 
& Carlberg 1991, Warren et al. 1992) , which produce 
non-axisymmetric disks if the halo dominates the inner 
mass density. Adding a dissipational gaseous component 
to these simulations leads to rounder, but still triaxial ha- 
los (Katz & Gunn 1991, Dubinski 1994). The intrinsic 
ellipticity and non-circular motions of these disks would 
produce large scatter in the TF relation (Franx & de Zeeuw 
1992) which is not observed. This limits either the halo tri- 
axiality to less than that predicted by any of the formation 
models, or the mass of dark matter within the optical disk. 
The shapes of the gravitational potentials of dark matter 
halos for two galaxies have been determined through dy- 
namical models of a polar ring (Sackett et al. 1994) and 
the flaring of H I gas (Oiling 1996). Both of these studies 
suggest that halos have very flat mass distributions, with 
axis ratios between 0.1 and 0.3, and thus resemble a disk 
more closely than a sphere. 

The various evidence discussed above strongly suggest 
that the luminous component contains a dynamically sig- 
nificant fraction (> 1/2) of the mass in the inner regions of 
spiral galaxies. However, most rotation curves (see Caser- 
tano & van Gorkom 1991) for two possible counter exam- 
ples) show no significant change in the transition between 
this inner region and the outer, dark matter dominated 
region. The disk and halo must somehow conspire to hide 
this transition (van Albada & Sancisi 1986). Moreover, 
because of the low scatter in the TF relation and in M/Ls, 
the disk-halo conspiracy must act consistently among all 
spiral galaxies to couple tightly the evolutionary history 
and the present structure of luminous disks and dark ha- 
los. This is extremely surprising, because after the initial 
collapse the evolution of the halo and disk proceed on dif- 
ferent scales and through different and complex physical 
processes which are linked only through a weak gravita- 
tional coupling. 

All of the above considerations led van Albada and San- 
cisi (1986) to advance the maximum disk hypothesis. Un- 
der this hypothesis, the mass of the luminous disk in a spi- 
ral galaxy is assumed to be as large as possible, consistent 
with the galaxy's rotation curve. The mass contribution 
of the dark matter halo is therefore assumed negligible in 
the inner parts of spirals. 

A maximal disk does not eliminate the disk-halo con- 
spiracy, in a sense, it makes it more puzzling because it 
minimizes the overlap in the distribution of dark and lumi- 
nous matter. However, it does make the tight correlation 
between mass and luminosity in the inner parts of spirals 
more plausible. 



Not all lines of evidence support the maximum disk hy- 
pothesis. Kuijken & Gilmore (1991) quote a local Galactic 
surface mass density, derived from the velocity dispersion 
of K dwarfs, which is 30% higher than that of "identified" 
matter. Cosmological N-body simulations of hierarchical 
universes yield halos which are not well approximated by 
isothermal spheres with cores (Dubinski & Carlberg 1991, 
Warren et al. 1992, Navarro et al. 1996). The mass den- 
sity profiles of these simulated halos have slowly changing 
logarithmic slopes and continue to rise all the way into the 
centers. The central mass concentration of these simulated 
halos would seem to be inconsistent with a maximum disk. 

In this work we model the luminous mass distribution 
for a sample of 74 spiral galaxies. Our aim is to test the 
maximum disk hypothesis by analyzing how well features 
in the rotation curve are reproduced by the mass mod- 
els. We therefore do not include a dark halo in the fits, 
but rather evaluate the quality of the fits under the strict 
maximum disk hypothesis. We extract surface bright- 
ness profiles for the disk and bulge from /-band images 
of normal spiral galaxies and derive optical Ha rotation 
curves from Fabry-Perot velocity maps with two spatial 
dimensions. The optical radius is defined as the radius to 
the extinction corrected 23.5 mag arcsec"^ isophote in /. 
The Fabry-Pcrot maps provide high signal-to-noise rota- 
tion curves, and enable us to average over local kinematic 
features. We fit axisymmetric mass models to optical rota- 
tion curves assuming constant, but distinct, mass-to- light 
ratios for the disk and bulge. The two dimensional infor- 
mation provided by the photometric images and velocity 
maps allows us to assess the importance of nonaxisymmet- 
ric features. This will be fully explored in a future paper. 

2. THE data/sample 

Our sample consists of 74 field and cluster spirals in the 
vicinity of the Hydra-Centaurus cluster. For each galaxy 
in the sample we have an /-band image and a two di- 
mensional Ha velocity map. Sixty-one galaxies are taken 
from Schommer et al. (1993) (SBWM). The observations 
and preliminary reductions of the data for these galaxies 
are described there. Observations and reductions for 13 
additional galaxies are presented in this paper. The sam- 
ple includes galaxies which are members of the clusters: 
Antha(6), Hydra(12), Centaurus(16), and Klemola 27 (5) 
(see Table 1). 

/-band CCD images were taken on March 15, 1994 with 
the CTIO 0.9m telescope at the Cassegrain focus. The 
detector was a TEK 1024 with a scale of 0.39"/pixel. The 
exposure time for each image was 10 minutes. Typical 
seeing was ~1.5". The images were bias subtracted and 
flatfielded with twilight skyflats using IRAF"^ . A large frac- 
tion of the SBWM images were taken with a TI 800x800 
CCD which had large, 1.5%, Hatfield errors. The errors 
were large scale and primarily near the edges of the field 
which made precise sky estimation difficult and created 
significant distortions at low surface brightness. An illu- 
mination correction was constructed for the TI images by 
combining all 40 /-band images from the run. For each im- 
age the galaxy and stars were removed and the remaining 
pixels were scaled by the sky level. The resultant image 
was smoothed over a scale of 20 pixels. The corrected 
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images are flat to 0.2%. 

All the images were calibrated using Graham standards 
(Graham 1982). The instrumental magnitudes of the 
standards were determined using DAOPHOT II (Stetson 
1987). DAOPHOT was also used to automatically locate 
cosmic rays and stars in the galaxy images. Stars and cos- 
mic rays were distinguished using the DAOPHOT sharp- 
ness parameter. The cosmic rays were removed by replac- 
ing affected pixels with the biweight of surrounding pixels 
(the biweight is an robust estimate of the mode of a dis- 
tribution which is optimized for small samples (see Beers 
et al. 1990). The stars were removed by flagging pixels in 
a circular aperture centered on each star. 

Ha spectroscopy for each galaxy was obtained in April 
1993 with the CTIO 1.5m telescope and the Rutgers 
Fabry-Perot imaging spectrophotometer. A TEK 512 
CCD with a 1.09"/pixel scale was used. The observations 
for each galaxy consist of 8 to 15 images with a 110 km s^^ 
FWHM bandpass (2.4 A at Ha), spaced at 1 A intervals to 
sample the Ha emission line. Wavelengths were calibrated 
during the day. Calibration drifts during the night were 
monitored by taking exposures of a neon lamp every hour. 
The maximum drift rate was 0.1 A/hr. The images were 
bias subtracted and flatfielded, with dome flats taken near 
the wavelength of each image, using IRAF. Typical seeing 
was ^1.5". 

Transparency and instrumental throughput variations 
were measured by performing photometry on stars com- 
mon to all the images for each galaxy using DAOPHOT. 
Cosmic rays were removed using the procedure described 
above. The images were convolved with a Gaussian to 
compensate for variations in seeing. The stellar positions 
were used to establish transformations and the images 
were shifted to a common coordinate system. The sky 
in each frame was determined in an annulus centered on 
the galaxy. There was generally no observable wavelength 
dependent structure in the the sky over the small 6 A 
gradient across Fabry-Perot field of view at the observed 
wavelengths. The series of images yields, at each pixel, a 
short segment of the spectrum around Ha. The spectra 
have been fitted with Voigt profiles (Humlicek 1979) to 
yield maps of the velocity, velocity dispersion. Ha inten- 
sity, continuum intensity, and their respective uncertain- 
ties. The kinematic data extends to a median of 4 disk 
scale lengths or 1.1 R23.5 in /. We used stars in each of 
the images to find relative astrometric corrections between 
the photometric and kinematic images. 

3. SPIRAL GALAXY MODELS 

We derive axisymmetric mass models of spiral galax- 
ies assuming that the radial mass distribution follows the 
radial luminosity distribution with constant M/L. The as- 
sumption of constant M/L requires that extinction and 
population gradients be small across the luminous disk, 
which may not be an adequate representation of real galax- 
ies, de Jong (1995) finds significant color gradients in the 
profiles of a large sample of spiral galaxies. His analy- 
sis attributes these gradients primarily to changes in stel- 
lar populations, with younger, more metal poor stars at 
large radii. This would tend to lower the stellar M/L at 
large radii and, therefore, decrease the radius at which the 
"missing mass" becomes important, de Jong's models pre- 



dict that (M/L)/ can change by a factor of 1.5 from the 
inner to outer parts of spirals. However, his models pre- 
dict a scatter of ~ 1.5 mag in the /-band TuUy-Fisher and 
a smaller scatter in color bands redder than /. Neither of 
these predictions are in agreement with observations. The 
large scatter is not observed and an analysis of the TF 
relation shows larger scatter both at visible color bands 
(Bothun & Mould 1987) and in the near infra-red (Bern- 
stein et al. 1994) than at /. These issues clearly need to 
be resolved, but in the absence of specific predictions for 
variable M/L we adopt a constant M/L for our models. 

We assume that luminous parts of spiral galaxies are 
composed of two principal components: a fiat disk and a 
rounder bulge. Each component is characterized by a dis- 
tinct spatial and kinematic stellar distribution, and stel- 
lar population. In addition, the disk harbors the cold gas. 
Our models are based on a two component, disk and bulge, 
photometric decomposition. To derive the mass distribu- 
tion we assume that each component has a separate, con- 
stant M/L. In real galaxies, the disk and bulge can be 
further divided into subcomponents such as the thin and 
thick disk, and nucleus; these are not considered here. The 
thick and thin disks in external galaxies have been distin- 
guished only in edge-on projections (Burstein 1979) and 
the nuclei are not resolved in our data. 

Spiral galaxies are approximately axisymmetric; the ax- 
isymmetry is broken by bars and spiral structure. Our 
models assume strict axisymmetry. The surface brightness 
profiles and the rotation curves for each galaxy are derived 
with a fixed and consistent center, position angle and incli- 
nation which we derive from both the photometry and the 
kinematic data. The difference in the projected disk and 
bulge axis ratios is taken into account in the disk-bulge de- 
composition. In this analysis there is no radial dependence 
of the geometric parameters. Such a dependence is often 
implicit in isophotal or kinematic tilted-ring analysis. 

Our goal is to examine how well maximum disk mod- 
els reproduce the mass distribution in a large and diverse 
sample of spiral galaxies. Individual models with halos 
could be constructed but would be poorly constrained by 
our data. 

In sections 3.1 and 3.2 we present our photometric and 
kinematic models. The 2-D nature of our data set allows us 
to derive independent geometric parameters from both the 
photometric images and the kinematic maps (SBWM). In 
section 3.3 we combine these results to derive a luminosity 
distribution and rotation curve with consistent geometric 
parameters. In section 4.3 we present the mass models. 

3.1. Photometric Models 

The goal of the photometric models is to separate the 
disk and bulge and deduce their radial luminosity distri- 
butions. The two primary features we use to distinguish 
the disk from the bulge in an image are the difference in 
the radial surface brightness profile and the axis ratio of 
the isophotes. 

The surface brightness profiles of most spiral galaxies are 
exponential over some fraction of the disk (Freeman 1970). 
Approximately 40% of spirals, however, deviate strongly 
from purely exponential disks. These spirals, designated 
as type II by Freeman, have profiles that are flat or slowly 
rising toward the center and have a steeper exponential 



outer profile. The profiles of type I galaxies more closely 
follow the canonical exponential disk. The profiles of disk 
galaxies often also exhibit smaller scale deviations. In or- 
der to model fully the radial structure of the disk these 
features must be included. We use an exponential com- 
ponent for the disk in our disk-bulge decomposition, but 
the final model for the disk is taken as the radially binned 
surface brightness profile of the galaxy after subtracting 
the bulge model. 

The bulges of spiral galaxies arc three-dimensional 
structures. This leads to the difference in axis ratios of 
the disk and bulge isophotes in the projected image of an 
inclined spiral. The isophotes of the bulge are typically 
not circular, and thus bulges arc typically not spherical. 

3.1.1. Disk-Bulge Decomposition 

There have been two strategies for performing disk- 
bulge decompositions. One has been to first perform 
an isophotal analysis; the resulting radial luminosity pro- 
file is fitted with an exponential disk and a bulge model 
(Schombert & Bothun 1987) such as the R^ law (de Vau- 
couleurs 1948) or Plummcr model. The second method, 
due to Kent (1986), uses the geometric properties of the 
disk and bulge. 

The profiles derived from isophotal fits are functions of 
up to five variables at each isophote: semi-major axis, el- 
lipticity, position angle and center coordinates. The de- 
pendence on the last four variables, while included in the 
fits, is often not presented. This dependence leads to dis- 
tortions in the radial luminosity profile which arc not in 
the underlying radial light distribution. The geometric 
parameters can be distorted locally by bright star forming 
regions. Unavoidably, the the parameters have a radial 
dependence that follows a bar and/or spiral arms. The 
isophotes are not well approximated by ellipses in the 
region where the surface brightness changes from bulge 
dominated to disk dominated. An elliptical isophote fit- 
ter makes a compromise fit which biases the true radial 
profile. For example, in a highly inclined galaxy with a 
compact nearly spherical bulge, the surface brightness of 
the bulge projects to large galactic radii along the minor 
axis. An elliptical isophote fitter will produce a profile 
in which the bulge seems to extend to larger radii. This 
will bias bulge models such as the R^ or Plummcr model 
which have strong tails. Byun & Freeman (1995) present 
a systematic study of these effects in model galaxies. 

The method of disk-bulge decomposition developed by 
Kent uses the different radial scaling properties of the disk 
and bulge profiles along the major and minor axes. The 
method relies on the different axis ratios of disk and bulge 
isophotes and the assumption of axial symmetry. The dis- 
advantage of this method is that it uses information only 
near the major and minor axis and is, therefore, sensitive 
to non-axisymmetric structure in these regions. 

Our method combines the best properties of these two 
methods by using the full two-dimensional information in 
the galaxy image. We fit two-components, a disk and a 
bulge model, to the image of the galaxy. We assume that 
the underlying distribution of light in spiral galaxies is 
axisymmetric, with a thin disk and an oblate spheroidal 
bulge with constant axis ratio. Under this assumption 
the projected isophotes of the model disk are ellipses with 



constant axis ratio, position angle and center. The pro- 
jected isophotes of the model bulge are ellipses with the 
same position angle and center as the disk but a differ- 
ent, larger, axis ratio. We exploit the exponential form 
of the disk but do not constrain ourselves to a more spe- 
cific functional form for the bulge. For the bulge we use a 
series expansion of Gaussians (Bendinelli 1991) which we 
have generalized to model oblate distributions. At large 
radii the surface brightness of this bulge model is always 
negligible compared to the disk. 

The disk and bulge are separated in three steps: (1) an 
exponential disk model is fitted to and subtracted from 
the image, (2) the bulge model is fitted to the resultant 
image, and (3) the bulge model is subtracted from the 
original image leaving an image of the disk. The exponen- 
tial disk of step 1 is fitted in a region well away from the 
bulge and extrapolated into the central bulge dominated 
region. The extrapolation is carried out in different ways 
for type I and type II disks (see below). The purpose of 
fitting the exponential disk model is to derive a global in- 
clination, position angle, and scale length and to provide a 
reasonable extrapolation of the disk into the central bulge 
dominated region for the purpose of isolating the bulge 
fight. 

For type I disks the galaxy image is divided into two 
regions; the disk dominated region and the central bulge 
dominated region. The disk dominated region is defined 
as an elliptical annulus centered on the galaxy: the inner 
edge of this annulus has axis ratio of the bulge isophotes 
and the outer edge has the axis ratio of the disk isophotes. 
The semi-major axis of the inner edge of the annulus is 
set sufficiently far from the center that the bulge is negli- 
gible compared to the exponential disk. The bulge dom- 
inated region is defined as an elliptical disk centered on 
the galaxy with the axis ratio of the bulge isophotes. The 
semi-major axis of each division is determined from a plot 
of the surface brightness binned in elliptical annuli with 
the approximate axis ratio of the disk. We fit a projected 
exponential disk to the image of the galaxy in the disk 
dominated region: 

l^o{a.^)^^^^e-''n/^n^ (1) 

where ii^ is central surface brightness, r^ is the disk scale 
length, and a^-, is the length of the semi-major axis of a 
ellipse of constant surface brightness for the disk. In gen- 
eral 

a^ = R'^{l+fsin'^{(l)-(f)o)) (2) 

where f^ — (-)^^ — 1 is defined in terms of the axis ratio, 
-, _R is the distance from the galaxy center as measured 
on the sky and 0o is the position angle of the major axis. 
For a flat disk the inclination is given by i = cos^^(-). 

The exponential disk is subtracted from the galaxy im- 
age and a bulge model is fitted in the bulge dominated 
region. The bulge model is a series of Gaussians: 
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where Ck is the total light in each component, r^ is the 

scale length of the fc*'* Gaussian component and a^ is the 
length of the semi-major axis of an ellipse of constant sur- 
face brightness for the bulge. Each Gaussian has a com- 
mon center and axis ratio. The position angle is fixed at 



the value derived for the disk, as required by the assump- 
tion of axisymmetry. Bulges are fitted with between 1 and 
6 Gaussian components. We use the maximum number 
of components that gives a unique and stable fit. If we 
add additional components we find either that: multiple 
components converge to the same scale length, with each 
of these components contributing a fraction of the inten- 
sity, or that the scale length of an additional component 
diverges and the intensity is reduced to the point that the 
component contributes a negligible constant offset. We 
have found that many of the standard fitting laws for the 
bulges of spiral galaxies such as de Vaucoulcurs or Plum- 
mer models are inadequate representations of real bulge 
images. The bulge model is subtracted from the original 
galaxy image, and the residual image is binned in ellip- 
tical annuli to arrive at the final disk model. To avoid 
truncation features and to integrate total magnitudes we 
extrapolate the disk model to large radii using the param- 
eters of the exponential disk. 

For type II galaxies we divide the galaxy into 3 regions: 
the outer disk region, the inner disk region and the central 
bulge region. We fit an exponential disk in the outer re- 
gion to derive the ellipticity and position angle. We then 
fit an exponential disk in the inner disk region with its own 
intensity and scale length, but with the ellipticity and po- 
sition angle fixed at the value derived from the outer disk. 
Because the surface brightness is nearly constant in this 
region, the shapes of the isophotes are not well defined 
and the light distribution contains no information about 
the position angle or inclination of the galaxy. We subtract 
the inner disk model from the image and the bulge model is 
fitted in the central region as described above. The bulge 
model is subtracted from the original galaxy image and 
the residual image is binned in elliptical annuli to arrive 
at the final disk model. The disk model is extrapolated 
to large radii using the parameter of the exponential disk 
derived in the outer region. 

The disk-bulge decomposition is not generally iterated, 
except to adjust the border of division between the disk 
and bulge regions. The disk model near the center is a rea- 
sonable, but arbitrary, extrapolation. Attempts to make 
small improvements in the solution through iteration de- 
pend on the details of the assumed disk extrapolation. The 
advantages of our method are that it does not start with 
built-in biases of the isophotal surface brightness profile, 
and it sets additional, geometric, constraints on the disk- 
bulge decomposition and therefore exploits the available 
2-D information. 

3.1.2. 3-D luminosity distribution 

The three dimensional luminosity distribution of the 
disk is trivially related to the projected distribution. As- 
suming no internal extinction p^{r) = /i^(r)cosi S{z), 
where r is the distance from the center of the galaxy and 
S{z) is the Dirac delta function. The surface brightness 
decreases by a factor of cosi in the dcprojection. 

Under our assumptions the two-dimensional elliptical 
bulge surface brightness distribution can be uniquely de- 
projected to the three-dimensional spheroidal luminosity 
distribution via Abel's integral equations (Stark 1977). 
A two-dimensional elliptical Gaussian distribution de- 
projects to a three-dimensional spheroidal Gaussian dis- 
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where a^ is the deprojccted semi-major axis of the 
spheroid. 
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follows directly from the galaxy inclination and the axis 
ratio of the bulge projected on the sky. 

3.1.3. Photometric Corrections 

The luminosity profiles are corrected for internal extinc- 
tion, Aint, and Galactic extinction, Ae^t- Corrections for 
internal extinction by dust generally give the total fraction 
of light absorbed within the galaxy. The surface bright- 
nesses of galaxies in this sample are derived by assuming 
the extinction is uniform over the luminous disk. The 
problem is, however, considerably more complex. The 
actual distribution of dust in galaxies is not well known 
and scattering may be as important as absorption in the 
I-band. Furthermore, multi-color photometry of spiral 
galaxies by de Jong (1995) suggests that radial color gra- 
dients may be due primarily to population gradients. 

We adopt Amt = — l.Olog(-) given by Giovanelli et al. 
(1994). Similar results are found by Han (1992), Bern- 
stein et al. (1994), and WiUick et al. (1995). Galactic 
extinction, Aext in the _B-band, is taken from Burstein & 
Heiles (1978). Reddening between B and / is assumed to 
be E(B-I)=0.45. The median value of Galactic extinction 
for this sample is 0.12 mag. 

3.2. Kinematic Models 

We derive rotation curves from two-dimensional Fabry- 
Perot Ha radial velocity fields. We assume that the Ha- 
emitting gas is in an axisymmetric rotating thin disk. In 
polar coordinates the model, projected on the sky, is given 

by 
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(7) 
where i is the inclination, (po is the position angle of the 
projected major axis, Vcirc{T) is the circular velocity pro- 
file and Vsys is the systemic velocity. The disk center is 
an implicit pair of parameters in the model. The term in 
brackets is equal to the cosine of azimuthal angle in the 
plane of the galaxy measured from the major axis. 

The parameters of the kinematic model are derived by 
fitting to the 2-D data in concentric elliptical annuli using 
a Levenburg-Marquardt x^ minimization technique (Press 
et al. 1992). The covariance matrix at the x^ minimum is 
used to estimate the errors in the parameters. The errors 
in the kinematic center are generally larger that those of 
the photometric center. The primary reason for this is that 



the kinematic center is poorly constrained along the minor 
axis and couples to the systemic velocity along the major 
axis. The center was therefore fixed by the centroid of the 
continuum distribution. In each annulus we fitted Vdrc, 4'o 
and i. The global kinematic position angle and inclination 
are the average, weighted by the estimated errors, of these 
parameters from each ellipsoid. The final rotation curve is 
extracted with all of the geometric parameters fixed. The 
rotation velocity is estimated independently on each side 
of the minor axis. 

3.3. Geometric Parameters 

The photometric and kinematic models yield indepen- 
dent estimates of the major axis position angle and the 
inclination. We merge these results and rederive the mod- 
els using consistent parameters. 

The position angle is generally better constrained by the 
kinematic models. A distinct line of nodes delineates the 
position angle in the velocity map, while the photometric 
position angle depends on the average distribution of lumi- 
nosity around an annulus. The photometric position angle 
is therefore more easily biased by global non-axisymmetric 
features such as spiral arms. 

The inclination, however is more poorly constrained by 
the kinematic model for many galaxies. For rising rota- 
tion curves the slope of the rotation curve and the inclina- 
tion are degenerate parameters. In the extreme case of a 
rigid rotator the degeneracy between the rotation profile 
and inclination is complete and there is no independent 
kinematic information on the inclination. SBWM found 
some galaxies with large deviations between the photo- 
metric and kinematic inclinations. The deviations are pri- 
marily for galaxies with rising rotation curves. 

We therefore determine the position angle from the kine- 
matic model and the inclination from the photometric 
model. The models are then iterated with the position 
angle fixed in the photometric model and the inclination 
fixed in the kinematic model. The final values of the geo- 
metric parameters are fixed and used consistently in both 
the photometric and kinematic models. The values are 
tabulated in Table 1. 



weighted average of the rotation curve points where the 
rotation curve becomes flat. For rotation curves which are 
still rising we define the width to be twice the maximum 
rotation speed. 

3.5. Mass Models 

We assume that the mass distribution of a spiral galaxy 
follows the de-projected luminosity distribution with con- 
stant M/Ls for each component. For the disk, we use 
a Fourier transform method for computing the rotation 
curve of a fiat axisymmetric mass distribution given by 
Kalnajs (1965): 
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where F is the Gamma function and j = \/— T. 

For a spheroidal mass distribution the rotation curve in 
the symmetry plane is given by Binney & Trcmaine (1987). 
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For a Gaussian distribution this integral reduces to a de- 
generate hypergeometric series in two variables. 

No dark halo is included. We also do not include a gas 
component. H I may contribute approximately 10% of the 
mass within the optical radius of late type spirals (Broeils 
& van Woerden 1994) and thus further lower the allowed 
mass of the luminous component. H I gas disks in spiral 
galaxies have longer scale lengths than the stellar disks and 
therefore contribute to the mass distribution primarily at 
larger radii where the rotation velocity due to the stellar 
disk begins to fall off. 



3.4. Galaxy parameters 

In Table 1 we list the photometric and kinematic param- 
eters for galaxies in the sample. After fixing the geometric 
parameters of the models, we derive values of the cen- 
tral surface brightness (So) and scale length (rj) for the 
disk. Dust extinction and projection effects are included 
in the models. We quote rj in the outer parts of the disk. 
For Freeman type II galaxies, therefore, the central sur- 
face brightness and scale length are not directly related 
quantities. The radius of the 23.5 mag arcsec"^ isophote 
(R23.5) is measured directly from the photometric profiles 
and, therefore, is also corrected for dust extinction (as- 
suming the dust acts uniformly over the disk) and projec- 
tion effects. The magnitudes for each galaxy are derived 
by integrating the disk profile and adding the total bulge 
luminosity. We quote total magnitudes integrated to in- 
finity. We calculate the ratio of the total luminosity bulge 
to that of the disk (B/D). 

We define the velocity width, used in the TuUy-Fisher 
relation, to be twice the rotation speed, measured by a 



4. RESULTS 

The model rotation curves of the bulge and disk are fit- 
ted to the kinematic data by adjusting the M/L of each 
component; the best fit is derived by minimizing x^- For 
most galaxies this results in a fit that nowhere significantly 
exceeds the data. For galaxies which would require a halo 
to get an acceptable fit we limit the radius of the fit so 
that the model rotation curve does not exceed the rota- 
tion curve data and evaluate the quality of the fit out to 
this radius. 

In Fig. 1 we present the models. The upper panel for 
each galaxy shows the face-on surface brightness profile. 
The contribution of the bulge is indicated by the dashed 
lines. The central surface brightness of the disk is marked 
with a diamond. The lower panel shows the rotation curve 
and the fitted mass model. The rotation speed for the side 
receding with respect to the center is marked with crosses 
and that for the the approaching side is marked with open 
circles. The rotation curves for the model bulge and disk 
are traced with dashed lines. The full model rotation curve 



is equal to the disk and bulge models summed in quadra- 
ture and is indicated by an unbroken line. In addition, 
the radii of prominent features, such as bars or rings, are 
marked with a vertical dot-dashed line. A scale bar in 
arcsec and the inclination are given to gauge the effect of 
seeing. 

4.1. Morphology of Surface Brightness Profiles 

In addition to the relative brightness of the disk and 
bulge, the surface brightness profiles of spirals can be clas- 
sified by the morphology of their disks. Two common as- 
sumptions are that spiral galaxies have an exponential disk 
and that they have a universal constant central surface 
brightness (Freeman's Law 1970). The discovery of low 
surface brightness galaxies (Bothun et al. 1991, Schombert 
et al. 1992, McGaugh et al. 1995) most clearly indicates 
that spiral galaxies fail the second assumption. The sur- 
face brightness profiles of disks are also not strictly expo- 
nential. The most prominent distinction in disk profiles is 
between Freeman types I and II, described in section 3.1 
Galaxy disks also frequently exhibit less extreme devia- 
tions from a pure exponential, such as a point of inflection 
where the scale length changes. 

One statement of Freeman's Law is that the ratio of 
the optical disk radius to the exponential scale length, 
R-23.5/rd, is a constant. In Fig. 2 we plot R23.5/rd vs. 
the central surface brightness. The tight correlation for 
exponential (i.e. type I) disks follows trivially from the 
definition of the optical radius and the the wide distribu- 
tion of central surface brightnesses in our sample. More 
significantly, R23.5/rd correlates with absolute magnitude 
(Fig. 3). We test the significance of this correlation with 
the Spearman rank-order correlation coefficient (Press et 
al. 1992). The correlation coefficient is —0.29. The prob- 
ability that there is no correlation is 0.01. The low lumi- 
nosity galaxies have lower central surface brightness and 
relatively flatter surface brightness profiles over the optical 
disk as indicated by a larger R23.5/rd. This leads to model 
disk rotation curves which rise more slowly for low lumi- 
nosity galaxies than those for high luminosity galaxies. 

4.2. Morphology of Rotation Curves 

Within the optical radius spiral galaxy rotation curves 
are not generally flat, but span a range of morphologies, 
from rising linearly to falling with radius (Rubin 1985, 
Persic & Salucci 1995). These authors also show that the 
shape of the rotation curves correlate with luminosity; low 
luminosity, small rotation velocity galaxies have rotation 
curves that are rising while high luminosity, large rotation 
velocity galaxies have falling rotation curves. 

An examination of the rotation curves in this sam- 
ple suggests such a trend, but a large range of rotation 
curve shapes is found at every scale of rotational veloc- 
ity. At high maximum rotation velocities (~300 km s^^): 
ESQ 374G02 (Fig. Iz) is falling, ESQ 375G12 (Fig. lac) 
is flat, ESQ 269G61 (Fig. li) rises over ~0.4R23.5 be- 
fore flattening, and ESO 381G51 (Fig. lai) rises over the 
entire optical radius. At low maximum rotation veloci- 
ties (~100 km s^i): ESO 374G03 (Fig. laa) is flat past 
~0.25R23.5, ESO 322G19 (Fig. Ik) is flat past -O.5R23.5, 
Abell 1644d83 (Fig. la) rises over the entire optical ra- 
dius. The rotation curve of ESO 441G21 (Fig. Iba) also 



rises to about 100 km s ^over the optical radius but does 
so almost linearly. 

4.3. Model Fits 

In spite of the variety of surface brightness proflle and 
rotation curve morphologies, maximum disk models with 
constant M/Ls provide good flts to optical rotation curves 
for a majority of the galaxies in our sample. The over- 
all structure of the rotation curves is reproduced by the 
models. The models fit the data for rotation curves that 
are rising linearly or with a curve, that are fiat or falling, 
or that have strong infiection points. The galaxies with 
good fits span a range of velocity widths (2wcirc) from 180 
km s-i for AbeU 1644d83 (Fig. la) to 680 km s^^ for 
ESO 572G17 (Fig. Ibv). For most of the galaxies no halo 
is required for the models to fit the data within the optical 
radius or to the last measured point. 

Fig. 4 shows the distribution of the ratio of the maxi- 
mum radius out to which the mass models provide a good 
fit to the optical radius; R/it/R23.5. The histogram is 
shaded for models which provide a good fit out to the last 
measured point of the rotation curve. 75% of galaxies for 
which the rotation curve is measured to R23.5 or beyond 
are well fit by a mass-traces-light model for the entire re- 
gion within R23.5. For 21% of the galaxies the models 
provide poor fits due to strong bars or spiral arms, these 
cases are assigned R/it/R23.5= in the histogram. The 
existence of strong non-axisymmetric structures suggests 
that there should be less dark matter in these galaxies. 
These cases are discussed below in section 4.4. 

For Freeman type I spiral galaxies, a thin exponential 
disk mass distribution has a rotation curve which reaches 
maximum at 2.15 disk scale lengths and falls slowly there- 
after, dropping 10% by 3.75 disk scale-lengths. Despite the 
apparent restriction of this shape, a wide range of optical 
rotation curve shapes can be successfully modeled because 
of the variation of the number of disk scale lengths within 
the optical radius among galaxies and the addition of a 
bulge component. 

Few galaxies have the canonical fiat rotation curve 
across the entire optical disk. ESO 375G12 (Fig. lac) and 
ESO 376G02 (Fig. lae) are two examples of galaxies with 
good fits that do. The optical radii in these galaxies span 
4.2 and 3.8 disk scale-lengths, respectively. The model 
rotation curve for ESO 375G12 begins to fall at the op- 
tical radius, but unfortunately, the rotation curve data 
for ESO 375G12 extend to only 3.3 disk scale lengths. 
ESO 376G02 shows slight evidence of dark matter, but 
only near the optical radius. The rotation curve for 
ESO 374G02 (Fig. Iz) has a gentle linear falloff over most 
of the optical radius and the model provides an excellent fit 
over this entire range. The optical radius of ESO 374G02 
extends over 4.1 disk scale lengths. These galaxies have 
prominent bulges. The flatness of the rotation curves is 
achieved though a combination of the bulge and disk ro- 
tation curves. 

ESO 383G88 (Fig. lam) and ESO 323g42 (Fig. Ix) have 
less prominent bulges. The rotation curves of these galax- 
ies rise less dramatically than those for the fiat rotation 
curves above, but do fiatten at larger radii. The op- 
tical radius of these spirals extend to 3.2 and 3.7 disk 
scale-lengths, respectively, and the models reproduce the 



turnover in the rotation curves. 

The optical disks of ESO 376G10 (Fig. laf) and ESQ 
501g01 (Fig. Ibl) extend to only 2.0 and 2.2 disk scale- 
lengths, respectively. These galaxies also have small 
bulges. The rotation curves reflect this morphology; they 
rise with a curve over most of the optical disk, reaching 
maximum near the optical radius. The rotation curve for 
ESO 501g01 extends significantly past the optical disk to 
4.1 scale lengths. The model rotation curve does not fall 
significantly below the data over this entire range. 

Freeman type II disks, which are distinguished by a flat 
inner surface brightness profile, are modeled with a con- 
stant inner mass density. This leads to rotation curves 
which are linearly rising over the constant density re- 
gion. The radius of the turnover varies from 20 to 60% 
of the optical radius for the galaxies in our sample. Ex- 
amples include: ESO 445G81, O.22R23 5 (Fig. Ibi); ESO 
375G02, O.23R23 5 (Fig- lae); ESO 216G20, O.3IR23 5 (Fig. 
Ic); ESO 501G86, O.55R23.5 (Fig. Ibp); and ESO 509G91, 
O.62R23.5 (Fig. lbs). The size of the linearly rising region 
in the rotation curve varies accordingly. 

ESO 381G51 (Fig. lai), ESO 435G26 (Fig. Ian), 
ESO 438G15 (Fig. law) and ESO 501gll (Fig. Ibm) also 
have flat surface brightness profiles near their centers; how- 
ever in these cases the profiles roll off slowly approaching 
an exponential asymptotically near the edge of the optical 
disk. The rotation curves in these cases rise slowly, with a 
curve matched by those of the models. ESO 435G26 has 
a strong bar and the fit within the bar radius is not good. 

The profiles of some spirals have inflections which are 
less prominent: the slope of the profile changes but does 
not become flat. This feature is reflected in the rotation 
curves of ESO 317G41 (Fig. Ij) and ESO 322G82 (Fig. 
Is). 

4.4. Nonaxisymmetric Structure 

A fraction, about 20%, of the fits fail in the inner re- 
gions; major structures in the model and/or rotation curve 
do not match up. These bad fits occur well within the op- 
tical disk and are not likely due to a dominant dark matter 
component. The galaxies which have the poorest fitting 
mass models often have strong non-axisymmetric struc- 
tures in the form of bars or strong spiral arm structure. 
These structures affect both the surface brightness pro- 
file and the rotation curve. The strong non-axisymmetric 
gas motions induced by a bar distort the measured rota- 
tion curve within the bar radius as noted for ESO 435G26. 
Strong spiral structure can affect the shape of the surface 
brightness profiles and induce large non-axisymmetric mo- 
tions as well as bias the determination of inclination and 
major-axis position angle. 

ESO 268g37 (Fig. Ig) and ESO 374g03 (Fig. laa) both 
have bars along the major axis and very strong spiral struc- 
ture. ESO 323g39 (Fig. Iw) has a bar along the minor 
axis and strong spiral ams. ESO 323g25 (Fig. lu) has 
very strong grand design spiral arms. The inclination of 
this galaxy, combined with the pitch angle and position 
angle of the spiral arms conspire in such a way that the 
arms closely follow the ellipticity of the disk from 5 to 10 
kpc. The large structure seen in the model at these radii 
is due to this chance alignment. This effect is also seen in 
ESO 322g36 (Fig. 11) and ESO 569G17 (Fig. Ibu) which 



also have strong spiral structure. 

A detailed analysis of the effect of strong non- 
axisymmetric structures will be presented in future paper. 
However, if a dark matter halo dominates the mass within 
the optical radius of these galaxies, it should act to stabi- 
lize the disk against these non-axisymmetric modes. It is 
difficult, therefore, to attribute the poor fit of these models 
to a large dark-to-luminous mass fraction. 

4.5. Small-Scale Structure 

Mass models for spiral galaxies sometimes, but not al- 
ways, reveal a correlation between small scale "bumps and 
wiggles" in the surface brightness profile and the rota- 
tion curve (Kent 1986, Freeman 1992). We find that in 
the two dimensional maps the residuals in the photome- 
try and the residual kinematic motions are highly corre- 
lated. How these structures show up in one dimensional, 
radial, surface brightness profiles and rotation curves de- 
pends strongly on how the data is sampled. The "bumps 
and wiggles" correlation is therefore most probably due 
to local perturbations, such as spiral arms and spiral arm 
streaming motions, rather than the global mass distribu- 
tion. 

4.6. Mass-to-Light ratios 

The median /-band M/L (disk plus bulge) for our sam- 
ple of galaxies is 2.4h75 in solar units with an rms scatter 
of 0.9. This accords well with the M/L predicted from stel- 
lar population synthesis models. If we exclude the galax- 
ies with the worst fits the median M/L rises to 2.7 ± 0.8. 
Worthey (1994) estimates that a normal stellar population 
which forms in a single burst will have an initial /-band 
M/L of about 1 in solar units and will increase to about 
5 over a time of 15 Gyr. His models have an M/L of 2.4 
at an age of 6 Gyr, and an M/L of 2.7 at 8 Gyr. Spiral 
galaxies are composed of stellar populations that span a 
range of ages and their M/L depends on the number of 
stars formed in each generation. A typical spiral might be 
about 10 Gyr old and form most of its stars over the first 
4 Gyr. 

The theoretical M/Ls are highly dependent on the as- 
sumed mass function for the stellar population and the 
details of the star formation history. Thus the agreement 
of our M/Ls with theory, although noteworthy, is not in it- 
self overwhelming evidence for the maximum disk hypoth- 
esis. The light comes primarily from a small population of 
bright, high mass stars, while the most of the mass is in a 
large population of faint, low mass stars. 

Current star formation rates vary with Hubble type, 
with lower rates in early type spirals. M/Ls increase over 
time as the stellar population fades, and therefore M/Ls 
should also vary with type, with higher M/Ls in early 
types. Such a correlation was found by Rubin (1985) and 
Kent (1986) for M/Ls in the F-band.' In the /-band this 
correlation is expected to be considerably weaker, and in- 
deed it is not strong in our data (Fig. 5). The stellar popu- 
lation models predict a rate of evolution of M/Ls in /-band 
that is a factor of two less than that in F-band for a given 
generation of stars. Also, the scatter in our M/Ls for a 
given type is large. 

The plot of M/L versus axial ratio (Fig. 6) is flat which 
gives us confidence that the assumed internal extinction 



law is correct and confirms the work of Han (1992), Gio- 
vanelli et al. (1994), Bernstein et al. (1994) and Willick 
et al. (1995). The good agreement of the values of M/L 
with the stellar population models also indicate that in- 
ternal extinction cannot be grossly higher, as claimed by 
Valentijn (1990). 

The scatter in the M/Ls in our sample is 37%. One 
of the primary sources of error is uncertainty in the dis- 
tance; measured M/Ls are inversely proportional to the 
assumed distance. The sample is concentrated in the great 
attractor region (Dressier et al. 1987), which is dominated 
by the Hydra-Centaurus superclustcr. The clusters are 
likely to have large peculiar velocities. Galaxies in the 
vicinity of the Centaurus cluster have a median M/L of 
1.8 ± .6 and galaxies near Hydra have a median M/L of 
3.4 ± 1. The Centaurus cluster is composed of two ma- 
jor sub-clusters (Lucey et al. 1986) and the Hydra cluster 
is also thought to have substructure (Fitchett & Merritt 
1988). The galaxies which are members of these clusters 
may therefore have especially large peculiar velocities. Fig. 
7 gives the distribution of M/Ls for the sample. Smooth 
Hubble flow is, unfortunately, not a precise approximation 
for estimating the distances to the galaxies in this sample. 
The large scale peculiar motions may contribute 25% or 
more to the error budget (Bothun et al. 1992, Mathewson 
et al. 1992, da Costa et al. 1996). This can also be seen 
in the TF relation for this sample. 

Figure 8 shows the TF relation for the galaxies in our 
sample assuming Hubble flow distances. We plot the to- 
tal magnitude with closed symbols for Freeman type I and 
open symbols for type II galaxies. We also indicate the 
R23.5 isophotal magnitude at the end of the line connected 
to each symbol. The length of the line indicates the ex- 
trapolation from the isophotal to total magnitude. The 
fainter galaxies clearly require a larger extrapolation; this 
is because the fainter galaxies have lower central surface 
brightness as noted above. 

If the slope of the TF relation is assumed to be near 6 
and we consider the total magnitudes (Mathewson 1992, 
Bernstein et al. 1994) then the scatter around the TF line 
is 0.46 mag. The zero point of the relation is fixed so that 
the average deviation from the TF line is zero (this as- 
sumes no bulk flow for the sample) . If the deviations from 
the TF line are attributed solely to peculiar motions, the 
error in the Hubble flow distances has a scatter of 24%. If 
the slope of the TF relation is assumed to be 10 and we 
consider isophotal magnitudes (SBWM, Peletier & Willner 
1993) then the scatter about the TF relation is 0.75 mag 
which implies deviations from the Hubble flow distances 
of 40%. 

5. CONCLUSIONS 

We find that spiral galaxy mass models which assume 
the maximum disk hypothesis yield good rotation curve 
fits within the optical radius, for a variety of spirals with 
distinct surface brightness and rotation curve morpholo- 
gies. 75% of galaxies for which the rotation curve is mea- 
sured to R23.5 or beyond are well fit by a mass-traces-light 
model for the entire region within R23.5. It is particu- 
larly striking that spirals with very different disk surface 
brightness profiles, generically distinguished by Freeman 
types I and II, are well modeled under the maximum disk 



hypothesis. 

Freeman type II galaxies constitute a significant frac- 
tion of spiral galaxies which fail to meet the canonical as- 
sumption that all spirals have an exponential disk. This 
has certainly contributed to the large discrepancy in disk 
scale lengths published by different authors for the same 
galaxies (Knapen & van der Kruit 1992). Type II galax- 
ies can be further characterized by the size of the inner 
fiat region and turnover rate. They span a large range of 
velocity widths and show no correlation with type. They 
generally break the correlation of rotation curve shape to 
absolute magnitude found by Rubin (1985) and Persic & 
Salucci (1995). 

Galaxies for which our models fail to give good fits in the 
inner regions, ^ 20% of the sample, generally have strong 
features, particularly bars but also strong spiral arms, 
which break the assumption of axisymmetry. Smaller scale 
deviations in the surface brightness profiles and rotation 
curves, "bumps and wiggles" , can also often be traced to 
nonaxisymmetric features. 

These results show that, within the optical regions of 
most spiral galaxies, the radial mass distribution is tightly 
coupled to the luminosity distribution. This is a much 
stronger constraint than that due to global correlations 
such as the TF relation. It implies that either the mass 
of dark matter must be small within the optical radius or 
that the distribution of dark matter must be precisely cou- 
pled to the distribution of luminous matter. A dark halo 
which is independent of, and unresponsive to the luminous 
disk cannot dominate the mass within the optical radius. 
A fraction of the derived luminous mass could be traded 
for halo mass, but the luminous mass traces the overall 
features of so many and various rotation curves that this 
fraction could not reasonably be too large. 

Persic and Salucci (1995) have constructed synthetic ro- 
tation curves by averaging over 500 optical rotation curves 
from Mathewson et al. (1992). They fit each synthetic ro- 
tation curve with model rotation curves for an isothermal 
halo and an exponential disk in which the scale length is 
fixed relative to the optical radius. They conclude that 
spirals galaxies with rotation velocities of 150 km s~^ are 
over 40% dark matter within the optical radius, and that 
galaxies with rotation velocities of 100 km s^^ are over 
75% dark matter. The assumptions which lead to this con- 
clusion are that all spirals galaxies have the same R23.5/rd 
and that all spirals of a given luminosity have the same ro- 
tation curve shape. Our results strongly suggest that both 
of these assumptions are too simplifying and call their re- 
sults into question. We find that R23.5/rd is larger in low 
rotation velocity galaxies. This leads to more slowly rising 
rotation curves. Relaxing only this assumption consider- 
ably weakens their dark matter constraints. 

Cosmological N-body simulations of hierarchical uni- 
verses suggest a universal halo profile for mass scales from 
ranging from dwarf galaxies to rich clusters of galaxies 
(Navarro et al. 1996). The halo profiles are centrally con- 
centrated and for spiral galaxies they dominate the mass 
distribution at all radii. In the cold-dark matter models 
of Navarro et al. , a galaxy with a maximum rotation ve- 
locity of 300 km s~^ is 72% dark matter within the optical 
radius, a galaxy with a maximum rotation velocity of 200 
km s^^ is 90% dark matter, and a galaxy with a maxi- 
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mum rotation velocity of 100 km s~^ is 96% dark matter. 
Our results show that optical rotation curves in real galax- 
ies exhibit a variety of shapes and that these shapes are 
well modeled by the luminous distribution of matter. It is 
difficult to reconcile our results with a universal rotation 
curve in which the central attraction primarily due to a 
dark matter halo. 

The success of the maximum disk hypothesis in model- 
ing the mass distribution in the inner parts of spiral galax- 
ies implies that either the mass of dark matter has to be 
small or that its projected distribution must follow pre- 
cisely that of the luminous matter out to nearly the optical 
radius. 
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Table 1 Notes 

Column 1- ESO-Uppsala Catalog number for galaxies in the sample. The codes for Figs. 1 are also given. 

Column 2- Cluster membership. Cen^Centaurus. 

Column 3- Distance in megaparsecs assuming Hubble flow with respect to the cosmic microwave background frame, and 
Ho= 75 km s"i Mpc~^ 

Column 4- Hubble type in RC3 (de Vaucouleurs 1991). 

Column 5- Freeman type. 

Column 6- Total /-band magnitude. 

Column 7- Disk central surface brightness in /. 

Column 8- Major axis radius of the / = 23.5 mag arcsee^^ isophote in kpc. 

Column 9- Exponential disk scale length. We quote the scale length of the outer disk for Freeman type II galaxies. 

Column 10- Bulge to Disk luminosity ratio. 

Column 11- Inclination. 

Column 12- Major axis position angle measured North through East. 

Column 13- log of the circular rotation velocity width in km s^^. 

Column 14 /-band M/L of the disk in solar units. 

Column 15- /-band M/L of the bulge in solar units. 

Column 16- /-band M/L of the disk plus bulge in solar units. 
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Fig. 1. — Mass models. Top panels: /-band surface brightness profile. The bulge profile is given by the dashed line. The diamond indicates 
the disk central surface brightness. Lower panels: Maximum disk fits to the rotation curves. The dashed lines indicate the bulge and disk 
contributions. See text for a full description. 
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Fig. 2. — Disk central surface brightness vs. R23.5/rd the sohd symbols are for exponential, Freeman type I disks and the open symbols are 
for Freeman Type II disks with a flat central disk profile. 
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Fig. 3. — Absolute magnitude vs. R23.5/i'd. The symbols are the same as in Fig. 2. The slope of the best fit least squares line is —0.37. 
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Fig. 4. — Histogram of R/it/R23.5- The shaded regions indicate galaxies that have good fits out to the last measured point. Cases where 
Rj^ij/R23.5= indicate galaxies which have poor fits due to bars or strong spiral arms. 
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Fig. 5. — M/L/ vs. Hubble type. The symbols are the same as in Fig. 2. The Hubble types, classified in the RC3 (de Vaucouleurs et al. 
1991), are coded; Sa=l, Sb=3, Sc=5, Sd=7, Untyped=10. 
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Fig. 7. — Histogram of M/L/. The hatched regions give the distribution of M/L in Centaurus and Hydra clusters. 
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